We present a novel GPU simulation algorithm based on projective dynamics (PD) with the penetration-free guarantee. This is achieved by incorporating barrier constraints into the PD framework. A faster GPU solver, named aggregated Jacobi, and a new CCD algorithm are also designed to further push the simulation performance on the GPU.